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We investigate quark potential by considering meson exchanges in the two flavor Nambu-Jona- 
Lasinio model at finite temperature and density. There are two kinds of oscillations in the chiral 
restoration phase, one is the Friedel oscillation due to the sharp quark Fermi surface at high density, 
and the other is the Yukawa oscillation driven by the complex meson poles at high temperature. 
The quark-meson plasma is strongly coupled in the temperature region 1 < T/T c < 3 with T c being 
the critical temperature of chiral phase transition. The maximum coupling in this region is located 
at the critical point. 
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I. INTRODUCTION 

It is widely accepted that there are two QCD phase 
transitions at finite temperature and density, one is the 
deconfinement and the other is the chiral symmetry 
restoration. The only possible way to realize the phase 
transitions in laboratories is through high energy nuclear 
collisions. From the experimentally observed anisotropic 
collective flow at RHIC which is characterized by the 
perfect fluid, the surviving resonant states of quarks and 
gluons above the critical temperature T c of the phase 
transitions calculated in lattice QCD and effective mod- 
els, and the significant difference between the lattice cal- 
culated thermodynamic functions and the corresponding 
Stefan-Boltzmann limits at extremely high temperature, 
the created new state of matter at RHIC is a strongly 
coupled quark-gluon plasma [l| . With the free energy of 
a pair of heavy quarks calculated in lattice QCD [|J , the 
heavy quark potential can be extracted. The depth of 
the potential well decreases with increasing temperature 
and the dissociation temperature of the bound state J/ 1 ^ 
is about 2.7 T c Q. 

From Yukawa's idea that nucleon-nucleon scattering 
is through boson exchanges, the static nucleon potential 
is a Fourier transformation of the boson propagator at 
momentum = (ko = 0,k), 



V(r) 



d 3 k 

(2V)3 



r C/(0,fc 2 ). 



(1) 



The Yukawa potential is widely discussed in nuclear 
matter and recently extended to quark matter, espe- 
cially, its two kinds of oscillations, the Friedel oscilla- 
tion 0, H, H, 0, H, Q induced by the sharp Fermi surface 
at high density and the Yukawa oscillation [f| [l(| re- 
sulted from the complex poles of the exchanged bosons, 
are frequently emphasized. 

To understand the quark behavior in strongly coupled 
quark matter with quarks, gluons and their bound states 
as constituents, we investigate in this paper the light 
quark potential through meson exchanges in a quark- 
meson plasma controlled by chiral dynamics. By calcu- 
lating the quark potential around the chiral phase transi- 
tion, we hope to see the effect of chiral symmetry restora- 



tion on the quark properties, to extract the information 
on the temperature region of strong coupling, and to 
know where the maximum coupling is. 

In our calculation we construct mesons in quark matter 
in the frame of Nambu-Jona-Lasinio 11] (NJL) model at 
quark level [l2j . The NJL model is inspired from the elec- 
tron superconductivity. It describes well the chiral prop- 
erties in vacuum and in hot and dense medium. Above 
the critical temperature, there exist hadronic excitations 
with small widths [HI, The spatial dependence of 

static meson correlation functions in the model was in- 
vestigated at finite temperature [l5| and density [l6l ]. 
In recent years, the model is successfully used to study 
color superconductivity at moderate baryon density [171 ] 
and pion superfluidity at finite isospin density [18j. To 
understand the results of lattice QCD thermodynamics 
in terms of quasiparticle degrees of freedom, the model 
is recently extended to include Polyakov loop dynamics 

(pnjl) Gi 13. 

To define the NJL model completely as an effective 
theory, a regularization scheme must be specified to deal 
with the important integrals that occur. A regulariza- 
tion scheme specifies a length scale for the theory, which 
is normally expressed as a cutoff A on the quark mo- 
mentum and meson momentum. From the uncertainty 
principle, the length scale in the NJL model is R ~ 1/A. 
Taking the standard cutoff value A ~ 600 MeV Q3, we 
have R ~ 1/3 fm. This means that, the quark potential 
calculated in the NJL model might be reasonable in the 
region of r > 1/3 fm. In the short range the model is 
probably not applicable. While we will take a covariant 
Pauli-Villars regularization scheme which allows us to do 
the momentum integrations in the whole region without 
an explicit cutoff, we still focus only on the long range 
behavior of the quark potential with r > 1/3 fm. 

The paper is organized as follows. We review the ana- 
lytic properties of meson polarization functions in the 
Pauli-Villars regularization scheme and determine the 
meson poles in the complex momentum plane in Sec- 
tion [il] In Section IIIII we show the quark potential at 
finite temperature and density and analyze the effect of 
chiral symmetry restoration. The temperature region of 
strongly coupled quark-meson plasma and the position 
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corresponding to the maximum coupling will be discussed 
in this section. We conclude in Section ITVl 



II. QUARK POTENTIAL 

The flavor SU(2) NJL model is defined through the 
lagrangian density 

C = ip (i^d^ - m + M7o) ip + G (i>^) 2 + (ip~ij 5 Tip) 

(2) 

with scalar and pseudoscalar interactions corresponding 
to a and 7r excitations respectively, where to is the cur- 
rent quark mass, G is the four-quark coupling constant 
with dimension (GeV)~ 2 , Tj (i — 1,2,3) are the Pauli 
matrices in flavor space, and [i is the quark chemical po- 
tential which is one third of the baryon chemical potential 
Ms- 

In mean field approximation the inverse quark propa- 
gator can be written as 



5- 1 (P)= 7 ^P M -M + 7 V, 



(3) 



where the dynamically generated quark mass M in the 
medium is related to the order parameter (ip^) of the 
chiral phase transition through M — m — 2G{tpip) and 
is determined by the self-consistent gap equation, 



M = to + 2iG 



d*P 



Tr [S{P)] 



(4) 



where and in the following the trace is taken in color, 
flavor and Dirac spaces. After the summation over Mat- 
subara fermion frequencies in the imaginary formalism of 
finite temperature field theory, the gap equation can be 
reexpressed as 



M(l - 2Gh) = mo. 



(5) 



Since the NJL model is non-renormalizable, one should 
take a regularization scheme to avoid the divergence in 
the momentum integration ([4]). By taking into account 
the covariant Pauli- Villars regularization approach, we 
can separate the function I\ in ([5]) into a vacuum and a 
matter part [l5j |. 

j rvac i rrnat 

h = h + L i 

N 



rvac 

h ~ 2ir 2 



fi N r°° n 2 

71 j=0 a=± J 

with the definition of E^ 1 (j>) = Ej (p) ± \i and Ej (p) — 



p 2 + A 2 , the Fermi-Dirac distribution function f(x) 
1 {e x l T + l) at finite temperature T, and the number of 



subtractions N. In the following numerical calculations 
we take the same parameter values as used in Ref. [l5j . 
too = 8.56 MeV, G = 0.75/2 fm 2 , N = 3, A = M, 
Ai = 680 MeV, A 2 = 2.1 Ai, A 3 = 2.1 A 2 and A = 1, and 
Ai, A2, A3 are determined by the constraints X) jLo A? 
fi. E^o^A 2 = and £^0^4 = 0. With these 
values one can fit the pion mass M w = 138 MeV and 
pion decay constant f n = 94 MeV and fix the quark 
mass M = 376 MeV and sigma mass M a = 760 MeV in 
the vacuum. 

In the NJL model, mesons are considered as quantum 
fluctuations above the mean field and can be constructed 
in random phase approximation (RPA) . After the bubble 
summation in RPA, the meson propagator is expressed 

as m 



U m {ko , k 



-2G 



i-2Gn m (fc 2 ,P) 

with the meson polarization functions 



(7) 



n m (fco > k ) — i 



d 4 P 



(2tt)4 

and the meson vcrtexes 



Tr 



r* m s(p + j)r n s(p 



2 ' 



(8) 
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m 


= a 




f 1 


m 


= a 




m 


= n+ 


r* = < 




m 


= 7T+ 




m 


= 7T_ 


m 




m 


= 7T_ 


*T375 


m 


= ttq , 






m 


= 7T0 



(9) 

Since the Lorentz invariance is explicitly broken at finite 
temperature, the polarization functions depend on fc 2 , and 
k 2 separately. 

With the NJL meson propagator (JT)) , the quark poten- 
tial through a meson exchange is expressed as 



V m (r) 



d 3 k 



(2tt) 3 
G 



e lk ' r C/ m (0,fc 2 ) 



2tt 2 ? 



Im / dk 



ke lkr 



1 - 2Gn m (0, fc 2 ) 



(10) 



By constructing a contour in the complex meson mo- 
mentum plane, the above integration can be done via 
calculating the meson pole in the complex plane. To this 
end, we need to know the analytic structure of the meson 
polarization function IL^fc 2 ,, fc 2 ). 

In the Pauli- Villars regularization scheme we again di- 
vide the polarization function into a vacuum and a matter 
part, 



n m (fc 2 , fc 2 ) — n™ c (fco, fc 2 ) + nj^ at (fc 2 , fc 2 ). 



(ii) 



From the discontinuity of the matter part on the real 
axis and the discontinuity of the vacuum part on the 
imaginary axis fl5j | . 
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Ren:: 



ReIC c (0,(ifc±e) 



rvac 



ReIC c (0,(fc± l£ ) 2 ) 

Tmat (0,(fc±ie) 2 ) = V{ 



3(fc 2 



2vr 2 



1 + fc? LaU/1 



(0 



fc 2 + % 



In A, 



mat _ 3(fc + £ m ) \ -t \ "» . 

2^ 2 fc ^ ^ j 



j=0 a=± 



dp 



E 3 {P) 



In 



k - 2p 



k + 2p 



f (£») . 



Imn™ c (0,(fc±*e) 2 ) = 0, 
Imnr t (0,(fc± l£ ) 2 ) = T 



3T(P + e 2 n ) 



N 



Akir 



i=0 a=± 



jvac 



2^ ^i=0 71 J 



N 



J- « / 1 — /c 2 arcsin /fy 



In A, 



Ren™ at (0,(i/c±e) 2 ) = /J 



ma I 



H4 



2ir 2 k 



EE^ 

j=0 a=± 



ImIC c (0,(^±e) 2 ) 
lmn^ at (0,(ifc±e) 2 ) = 



± 



3(e 2 m -k 2 ) 
iir 



k) 



dp 



1 In fc 



In A,, 



for kj < 1 
for kj > 1 



Ei(p) 



7r — 2 arctan 



for fcj < 1 
for fcj > 1 



(12) 



r 



with the definition of kj = k/(2Aj) and e 2 ^ = for m = 
7T + , 7r_, 7To and e 2 n = AM 2 for m — a, there is a vacuum 
cut on the imaginary axis starting at 2m and a matter cut 
on the real axis in the complex meson momentum plane. 
To replace the direct momentum integration p0[) by the 
contour integration in the complex momentum plane, we 
need to know the pole position of the meson propagator. 
Suppose the pole is at k = iM m + T m , the mass M m and 
width T m are determined by the complex pole equation 



1 - 2Gn m 0, {iM m + r„ 



0. 



(13) 



Since the Lorentz invariance is broken at finite temper- 
ature, the mass M m and width T m related to the quark 
potential are in general different from the dynamic meson 
mass M m and the corresponding width T m controlled by 
the pole equation 



i-2Gn m ((A/ m -zr m ) 2 ,o) = 0. 



(14) 



Only in the vacuum with T = jl = 0, the two poles 
coincide. 

Suppose the width is much less than the mass, T m << 
M m , the mass is decoupled from r m , 



1 - 2G Ren m (0, (iM m + ef) = 0, 

and r m is determined by 

e 2 m -Mj Imn m (0,(»M m + e) 2 ) 
2M m Ren m (0, (*M m + e) 2 ) - h ' 



(15) 



(16) 



Due to the mass constraint M ff < 2M at T < T c , pions 
are bound states and can not decay into two quarks, and 



therefore they have no width in the chiral breaking phase. 
In chiral limit with zero current quark mass, from the 
mass relations M m — e m at T < T c , not only pions but 
also sigma have no widths in the chiral breaking phase. 

Taking into account the contributions from the vacuum 
cut, matter cut and the pole to the total quark potential, 
we have 

V m (r) = VZ ac (r) + Vr\r) + V^ l {r) (17) 
with the vacuum part 



VZ ac (r) 



G 
2n 2 r 



dkke Im 



2 M 



1 



l-2GIT m (0, (ik + e) 



1 



l-2Gn m (0,(ifc-e) 2 ) 
which decays exponentially, the matter part 
G 



(18) 



vr\r) 



2n 2 r 



dkk Im e 



1 



■!•■)■ 



1 



l-2GII m (0,fc 2 ) 
(19) 



1 -2GIi m {0,{k + it) 2 ) 

which in general oscillates, and the pole part 

V pol ( . _ 1 e ~ M ^r am cos (r m r) + b m sin(r m r) 
m [r '~2n r a 2 n + b 2 m 

with the coefficients a m and b m defined by 

dU m (0,k 2 ) 



dk 2 



(20) 



(21) 
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where the mass and width control, respectively, the am- 
plitude and the oscillation frequency of the pole contri- 
bution. In the chiral breaking phase, the pole is exactly 
on the imaginary axis, we recover the traditional Yukawa 
form 



0.8 



1 e 



-M„ 



2tt a r . 



(22) 



III. RESULTS 



We first show the effect of chiral symmetry restora- 
tion on the quark potential in hot and dense medium. 
With the parameter values in the model and in the Pauli- 
Villars regularization scheme, shown in the last section, 
the pole mass M m and width T rn as functions of temper- 
ature at zero chemical potential are presented in Fig[T] 
for m = 7r, a. In chiral limit, the critical temperature 
T c for chiral phase transition is about 300 MeV. With 
different parameter values or choosing other regulariza- 
tion schemes, the critical temperature can be reduced. 
For instance, in the three-momentum non-covariant cut- 
off scheme, the critical value is T c = 170 MeV p3 |. In the 
real world with nonzero current quark mass, pions and 
sigma behave differently in the chiral breaking phase. Pi- 
ons are bound states with zero width, but sigma mass is 
slightly larger than two times the quark mass. Therefore, 
there is a small width with sigma at low temperature. 
When chiral symmetry is restored at high temperature 
T > T c , the difference between pions and sigma disap- 
pears, both the masses and widths coincide. The width is 
always much smaller than the corresponding mass. Even 
at high temperature T/T c = 3 the ratio T m /M m is only 
about 10%. This is in consistent with the assumption of 
r m << M m for the decoupling of M m from T m , (fT3|> and 
(I16p . The meson masses get saturated when the temper- 
ature is high enough, this will lead to the saturation of 
the quark potential, see the following discussion. Note 
that, the finite width plays a crucial rule in the chiral 
restoration phase. At T > T c , quarks become almost 
massless, and pions and sigma satisfy the decay condi- 
tion M m > 2M. If the meson width is neglected, the 
pole is located on the imaginary axis and goes into the 
vacuum cut, and the pole contribution will be fully sup- 
pressed. However, for a nonzero width, even if it is very 
small compared with the mass, the pole position deviates 
from the imaginary axis, it is outside the vacuum cut and 
its contribution is not removed from the quark potential. 




FIG. 1: The quark mass, meson masses (upper panel) and 
meson widths (lower panel) as functions of temperature at 
zero chemical potential. The temperature is scaled by the 
critical value T c . 
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FIG. 2: The pion-mediated quark potentials with and with- 
out considering chiral symmetry restoration. The results are 
compared with the potential in the vacuum. 



To directly see the effect of chiral symmetry restora- 
tion on the quark potential, we show in Figj2] the pion- 
mediated potentials at T/T c = 1.3 and /j, = in the 
cases of with and without chiral symmetry restoration, 
and compare them with the potential in the vacuum 
(T = (J, = 0). For the calculation without chiral phase 
transition, we keep the quark mass and pion mass as 
their values in the vacuum and remain only the explicit 
temperature dependence in the distribution function /. 



In this case, the quark potential becomes stronger in 
comparison with the vacuum potential, because of the 
increasing quark density in hot medium, and is signifi- 
cantly different from the calculation with chiral symme- 
try restoration where the potential is much weaker than 
the vacuum potential, due to the change in the chiral 
dynamics. 

Since sigma is much heavier than pion in the phase 
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with chiral symmetry breaking, and they become degen- 
erate when the symmetry is restored, sigma exchange is 
remarkably reflected in the quark potential only in the 
short range for T < T c and becomes equivalently im- 
portant as pion exchange for T > T c . This is explicitly 
shown in FigO 
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FIG. 3: The pion- and sigma-mediated quark potentials in 
the chiral breaking and chiral restoration phases. 




FIG. 4: The pion-mediated quark potential and its three com- 
ponents arising from the vacuum cut, matter cut, and pole 
in the chiral restoration phase at high temperature (fi = 
and T/T c — 1.2, upper panel) and high chemical potential 
(T = 130 MeV and fi/fi c = 1.1, lower panel). 



We now discuss the two possible oscillations in the 
quark potential, the Friedel oscillation and the Yukawa 



oscillation. The former is due to the sharp Fermi sur- 
face and the latter arises from complex poles of meson 
exchanges. From the matter and pole contributions (TT51) 
and ptJ|) , there are always oscillations in the two parts. 
To know whether they can be clearly reflected in the to- 
tal potential, we separately show the three contributions 
from the vacuum cut, matter cut and pole in Fig[U While 
there is a remarkable oscillation around r = 1 fm in the 
matter part at fi = and T/T c = 1.3, see the upper 
panel, it can not be seen clearly in the total potential, 
since at high temperature and low density the matter 
contribution is not the dominant one in the total poten- 
tial. In this case, only Yukawa oscillation can be seen, if 
the width is large enough. At T/T c = 1.3 the pion width 
is too small and the induced weak oscillation is difficult 
to be demonstrated in the total potential. On the other 
hand, at low temperature and high density, the quark po- 
tential is absolutely controlled by the matter cut, see the 
lower panel at /Lt//x c = 1.1 with fi c = 400 MeV being the 
critical chemical potential for T = 130 MeV, and only 
the Fermi-surface induced Friedel oscillation can be seen 
in the total potential. When the temperature and chem- 
ical potential are not low enough, especially in the chiral 
restoration phase, the contribution from the vacuum cut 
is always smaller than the matter part and pole part. 

The Friedel oscillation is related to the singularity of 
the matter part of the meson polarization function II m 
at k = 2pf on the real axis at zero temperature, where 
Pf is the quark Fermi momentum, see (|12p . At finite 
temperature, there is no more strict Fermi surface and 
singularity in the polarization function, the oscillation 
will become weak at low temperature and finally disap- 
pear at high enough temperature. In the upper panel of 
Fig|5]we plot the quark potential at fixed chemical poten- 
tial for three values of temperature. The Friedel oscilla- 
tion is gradually washed away by increasing temperature. 
The chemical potential dependence of the Friedel oscil- 
lation at fixed temperature is shown in the lower panel. 
In the chiral breaking phase with low chemical poten- 
tial /i < M, there is no Fermi momentum and in turn 
no Friedel oscillation. The oscillation happens only in 
the chiral restoration phase, and the amplitude and fre- 
quency increase with increasing chemical potential. 

The strength of the Yukawa oscillation depends 
strongly on the pole width. In FigJB] we show the pion- 
mediated quark potential in the chiral restoration phase 
with [i — and T/T c = 2 for three different values of 
width. As we emphasized above, the potential at high 
temperature and low density is dominated by the pole 
contribution. From the analytic structure of the meson 
polarization function (|12p . its matter part is continuous 
on the imaginary axis, and only the vacuum part con- 
tributes to the width (|16[) . From our numerical calcula- 
tion, see Fig[Tl the width is only about 10% of the mass at 
T/T c = 2, and therefore the Yukawa oscillation can not 
be clearly shown in FigO However, this small width can 
not be neglected, otherwise the pole will be located on 
the imaginary axis and go into the vacuum cut, and has 
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FIG. 5: The temperature (upper panel) and chemical poten- 
tial (lower panel) dependence of the Friedel oscillation in the 
pion-mediated quark potential. 
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FIG. 6: The width dependence of the Yukawa oscillation in 
the pion-mediated quark potential in the chiral restoration 
phase. 



no contribution to the potential. To see a clear Yukawa 
oscillation, we take artificially the width as a free param- 
eter and calculate the potential with T^/M^ = 2. In this 
case, there is really a strong oscillation around r = 1 fm, 
as we expected. 

From the analysis of the experimental data obtained in 
relativistic heavy ion collisions, especially the anisotropic 
flow, and the lattice and model calculation on the reso- 
nant states of quarks and gluons above the critical tem- 



perature, it is widely accepted that the parton system 
above and close to the critical temperature is a strongly 
coupled quark-gluon plasma. The temperature region is 
about 1 < T/T c < 3 ~ 4 [2l[. Above this region, quarks 
and gluons are weakly coupled and the system can ap- 
proximately be considered as an ideal gas. In Fig[7| we 
show the temperature dependence of the pion-mediated 
quark potential at zero chemical potential. With in- 
creasing temperature, the potential becomes weaker and 
weaker due to the chiral symmetry restoration, and fi- 
nally gets saturated for T/T c > 3. Since the used Pauli- 
Villars regularization scheme guarantees the ideal gas 
limit at T — > oo (l5j . we can conclude from this satura- 
tion that, the quark- meson plasma is strongly coupled in 
the temperature region 1 < T/T c < 3. Above T/T c ~ 3, 
the system behaves like an ideal gas. While there are 
no gluons and confinement in the NJL model, the quark- 
meson plasma has the similar strong coupling region as 
the quark-gluon plasma. 




FIG. 7: The temperature dependence of the pion-mediated 
quark potential at zero chemical potential. The potential gets 
saturated at T/T c = 3. 



To describe the degree of particle coupling in matter, 
one can introduce the ratio of particle potential to parti- 
cle kinetic energy (22[, 



7 = V /E ki 



(23) 



A large ratio means a strongly coupled matter and a small 
ratio indicates a weakly coupled matter. For the consid- 
ered quark-meson plasma, we take the potential Vq as the 
meson-mediated quark potential at the average distance 
between two quarks, 



V = \V(2r )\, 



(24) 



where V is the summation of V m over all possible meson 
exchanges, the averaged quark radius ro is determined by 
the quark and anti-quark number densities n q and rig, 



(n q +n q ) -nr = 1, 



n q ,q(T, n) 



d 3 p 
(2VP 



f(E(p) T »), (25) 
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with E(p) = -\/p 2 + M 2 , and the quark kinetic energy is 

x (f(E(p) - /<) + /(E(f>) + rt) . (26) 

The ratio is shown in Fig[8] as a function of chemical 
potential at fixed temperature T — 130 MeV. Different 
from the kinetic energy E^m which monotonously goes 
up with increasing temperature or chemical potential, the 
potential Vo depends on both the shape of the potential 
and the quark density. The former is controlled by chiral 
dynamics and the latter is determined by thermodynam- 
ics. With increasing temperature or chemical potential, 
the quark potential becomes weak, while the density in- 
creases and in turn the radius ro gets small. From our 
numerical calculation, the ratio is rather small at low 
and high chemical potentials, but peaks strongly at the 
critical point. This tells us that the strongest coupling 
in the quark-meson plasma is located at the chiral phase 
transition. The smoothly increasing ratio after the phase 
transition is due to the strong Friedel oscillation at high 
density. From the AdS/CFT estimation [23| and pertur- 
bative QCD calculation [24|, it is well-known that the 
ratio of shear viscosity to entropy rj/s has the minimum 
value at the critical point of any phase transition, which 
is widely considered as a general property of strongly cou- 
pled matter. The maximum of the ratio Va/Ekin at the 
chiral phase transition agrees well with the minimum of 
the ratio tj/s, both describe the strongest coupling at the 
phase transition. 




0.20.40.60.8 1 1.21.4 



FIG. 8: The ratio of quark potential to quark kinetic energy 
as a function of chemical potential at T = 130 MeV. The peak 
is located at the critical point of chiral symmetry restoration. 



IV. CONCLUSIONS 

In the frame of the two flavor SU(2) NJL model, we cal- 
culated the light quark potential through pion and sigma 



exchanges in RPA approximation. The chiral symme- 
try restoration at high temperature and density plays an 
essential rule in the quark potential. While the mat- 
ter part is always associated with an oscillation, it dom- 
inates the potential only for dense matter, and there- 
fore the oscillation can be reflected clearly in the total 
quark potential only at low temperate and high density. 
The other possible oscillation in the chiral restoration 
phase is the Yukawa oscillation induced by the complex 
poles of mesons. From our self-consistent calculation for 
meson mass and width, the Yukawa oscillation is weak 
even at extremely high temperature. The potential gets 
saturated when temperature is larger than about three 
times the critical temperature. This means that the 
strongly coupled quark-meson plasma is in the region of 
1 < T/T c < 3, which agrees well with the widely accepted 
result for the strongly coupled quark-gluon plasma. Cor- 
responding to the well-known minimum of the ratio of 
shear viscosity to entropy, we found that the ratio of po- 
tential to kinetic energy peaks strongly at the critical 
point of chiral phase transition. 



While we found the strong coupling region for the 
quark-meson plasma through considering the saturation 
of quark potential, the saturation temperature T/T c = 3 
looks beyond the effective range of the NJL model. The 
reliable conclusion is probably that the quark-meson 
plasma above and close to the chiral critical point is 
strongly coupled. In our calculation we did not include 
the vector meson exchange which is expected to con- 
tribute a repulsive part to the quark potential and then 
reduce the depth of the potential well in the short range. 
To have a complete calculation of the quark potential in 
the NJL model, one should consider the vector mesons. 
Another point is the dcconfincmcnt effect, it will be 
helpful to investigate the quark potential using the im- 
proved NJL model with Polyakov loop dynamics [l^, H(| • 
Since the original NJL model with pions and sigma 
describes well the chiral dynamics in the vacuum and at 
finite temperature and density, we hope our result on 
the quark potential in the quark-meson plasma is helpful 
to understand the effect of chiral symmetry restoration 
on the strongly coupled quark-gluon plasma. 
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